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В статье анализируется влияние инструментальных погрешностей и погрешности квантования угловых 
перемещений на точность аналого-цифрового преобразования угловой скорости подсчетом числа единич- 
ных приращений кода углового перемещения за фиксированный интервал времени и определены 
точностные требования к выбору основных параметров цифровых тахометров, построенных на основе 
растровых синусно-косинусных преобразователей, и синтезированы точностные и динамические модели, 
позволяющие исследовать их в составе цифровых электроприводов. 
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Введение 


В безредукторных цифровых электроприводах прецизионных мехатронных ин- 
теллектуальных систем широкое применение получили датчики местной, скоростной и 
позиционной обратных связей на основе фотоэлектрических цифровых преобразова- 
телей перемещений, обеспечивающих достижение быстродействия и точности в рас- 
ширенном диапазоне изменения скорости вращения, минимальных массогабаритных 
размерах информационно-силового модуля и в условиях значительного изменения 
момента нагрузки [1], [2]. 

Из них наибольшее применение получили инкрементные преобразователи 
перемещений, построенные на растровых синусно-косинусных преобразователях пере- 
мещений (СКП), что объясняется их широким распространением и возможностью 
формирования с их помощью как позиционного, так и скоростного каналов преобра- 
зования в виде цифровых тахометров и тахометров постоянного тока, основанных на 
дифференцировании угла по времени в процессе аналого-цифрового, цифрового и 
аналогового преобразований [3-5]. 

Для регулирования скорости вращения на высоких и средних частотах враще- 
ния в основном используются цифровые тахометры, основанные на преобразовании 
частоты следования импульсов инкрементных преобразователей в цифровой код. Их 
достоинством является простота построения, а недостатком — низкая разрешающая спо- 
собность, зависящая от коэффициента умножения пространственной частоты растрового 
СКП, определяющего число формируемых импульсов за один оборот входного вала 
и соответственно коэффициент масштабного преобразования угловой скорости вра- 
щения [2-4]. 

С целью увеличения числа импульсов на оборот и расширения динамического 
диапазона регулирования скорости вращения до 10° и выше используются инкрементные 
преобразователи повышенной разрешающей способности в виде структур, состоящих 
из последовательно соединенных СКП, блоков электронного умножителя простран- 
ственной частоты (УПЧ) и формирователя импульсов единичных приращений (ФЕП) [4]. 

Однако вопросы их точностного анализа и синтеза недостаточно разработаны, 
что затрудняет выбор их оптимальной структуры при проектировании. 

Целью работы является анализ погрешностей, порождаемых инструменталь- 
ными погрешностями растровых СКП, погрешностями квантования угла и угловой 
скорости, и точностной синтез параметров цифровых тахометров. 


Анализ влияния инструментальной погрешности 
на точность преобразования угловой скорости 


Параметры скоростного канала определяются инструментальной погрешностью 
преобразователя угла Да», шагом квантования угла и скорости 0, и ©, а также 
шагом дискретизации отсчетов Та. Рассмотрим требования к выбору перечисленных 
параметров для случая дифференцирования угла по времени методом прямой 
функции преобразования (рис. 1), обеспечивающим постоянство шага дискретизации 
отсчетов во всем диапазоне преобразуемых скоростей © [6]. 

Схема позволяет осуществить преобразование скорости в код с разрешающей 
способностью [5] 


2п 
= 
Мы М 
= коэффициенты умножения пространственной частоты блоков 


> 


где М ,М№ ‚№ 


скп’ °'упч’ `` феп 
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СКП, УПЧ и ФЕП; Т, — образцовый измерительный временной интервал; № — 
преобразуемая частота импульсов инкрементного преобразователя перемещений. 


Рисунок 1 — Структурная схема цифрового преобразователя угла и скорости 


Преобразование угловой скорости в код данным методом осуществляется 
нахождением приращений кода углового перемещения за фиксированный интервал 


времени Т <Т,, равный образцовому измерительному временному интервалу, либо 
подсчетом импульсов единичных приращений перемещения за время интервала Т, либо 
вычислением разности кодов угловых положений в моменты его начала и окончания. 
При наличии Ло» погрешность преобразования скорости АО, может быть пред- 
ставлена с точностью до членов первого порядка малости суммой двух слагаемых 


АО, =ЛО, АО, (1) 
ах А, (№, )- А, (№0, ) | 8 
т 
АО, =-ЛЁ „©, (3) 
где ДО„ — инструментальная погрешность преобразования 62; & а {® а га -: 


шаг квантования © и <; Дб), — погрешность квантования (2; ЛЕ =&,-&,, М№,, — разность 


погрешностей квантования угловых положений, выраженных в Долях 0 ‚ И значе- 


ния кода угла, соответствующие моментам начала и окончания интервала Т. 
Подставив в (2) вместо До его разложение в ряд Фурье 


А, -Ул зш (Ки + у, ), 
К=1 


где о = ОЕ, будем иметь 


АО, = У За Е КОЛ, с0$| К| ©, ВЯ А (4) 
к 2 2 
зш (х) : 
где За(х) — функция ———^; Л,, у, — амплитуда и фаза А-й гармоники по- 
ь 


грешности Ао; г. — момент окончания временного интервала Т. 


й 

Анализ (4) показывает, что АФ„ представляет собой дискретизированную по 
времени с шагом Т производную инструментальной погрешности Ло; по времени, 
ОГ 


амплитуды гармоник которой умножены на ы } а фазы задержаны на 
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КОТ ь 
величину —, что соответствует прохождению инструментальнои погрешности 
2 2 


через дифференцирующий элемент и элемент с передаточной функцией [6] 


а 


Й (5) 


Принимая во внимание, что восстанавливающий элемент при ступенчатой 
аппроксимации отсчетов имеет передаточную функцию [6] 


а" 


Й =—_, 6 
(р) г (6) 


математическая модель скоростного канала может быть представлена в виде импульс- 
ной линейной системы (рис. 1 а), состоящей из последовательно соединенных диффе- 
ренцирующего, суммирующего элементов, элемента с передаточной функцией (5), клю- 
чевого элемента и элемента с передаточной функцией (6). 

Импульсная система (рис. | а) может быть заменена непрерывной (рис.1, 6), если 
частота дискретизации превышает удвоенную максимальную частоту Лош. В этом 
случае импульсный элемент может быть представлен пропорциональным элементом 


с коэффициентом передачи г 


С учетом этого погрешность преобразования скорости ДО, может быть пред- 
ставлена выражением 


А=К тах 
АО, = У, ба’ [Г вол, соз[ © (#-Т)+чи, |. (7) 


К=1 
Анализ (7) показывает, что наиболее сильное влияние Ло„ будет сказываться 


к 


на низких скоростях, где м =1. При этом относительная погрешность будет равна 


—_Ющах 


КК, 
0, = У, КА, соз[ «0 (1-Т)+ чу, |. (8) 
К-1 


Аои(р) 


6) 


Рисунок 2 — Математическая модель скоростного канала 
цифрового преобразователя угла и скорости без учета погрешности квантования: 
а) представление скоростного канала линейной импульсной системой; 
6) представление скоростного канала линейной непрерывной системой 
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Выражение (8) позволяет оценить требования к характеру погрешности пре- 
образователя угла, исходя из допустимой погрешности ДО, на низких скоростях <> 


осей < м : (9) 
фе {® 


Анализ влияния квантования угла на точность 
преобразования угловой скорости 


Представив (3) в следующем виде 
_ Аа, (Г) - Аа, (АТ-Т) ий 
«_ Т 2 

нетрудно видеть, что погрешность квантования угловой скорости А, обусловлена 
прохождением через звенья импульсной схемы (рис. 2, а) погрешности квантования 
угла Ло.. Поэтому математическая модель скоростного канала с учетом погрешности 
квантования угла может быть представлена импульсной линейной системой, изобра- 
женной на рис.З3. 


[о 
Принимая во внимание, что Да, пена а = ‚где Е т| ] — дробная часть числа, 
Ч 


заключенного в квадратные скобки, из (10) находим 


АО, =-0, ть | (10 


где Ао -— погрешность квантования в начальный момент времени. 
Представив зависимость угла во времени в виде © (1 ) = О+ Ао, (© ) и разложив 


правую часть (11) в ряд Фурье, будем иметь 


ло = ей У я, (о/+ а, (9/)+А,)- 


Ч 
ды 


тим, (©(!-Т)+ Аа, (©(и-Т))+ 4, |. (12) 


2п я 
где №, = —— - количество уровней квантования угла в скоростном канале. 


Рисунок 3 - Математическая модель скоростного канала цифрового 
преобразователя угла и скорости с учетом погрешности квантования: 


Ла, (р) — изображение функции До, (1) по Лапласу 


416 «Искусственный интеллект» 2013 №3 


Анализ погрешностей преобразования и точностной синтез... 


Свернув разность синусов в квадратных скобках выражения (12), получим для 
низких частот ©, на которых влиянием элемента с передаточной функцией (5) можно 
пренебречь, 

ь гу, (1+ да, (01) т 
49, =29(1+ Аа, (©) У За г 
1=1 


о(!+ ла, (01) т 
хсо$ У, фича, (и Е (13) 


р алое, ы 
где Ла, (©) = о производная инструментальной погрешности по углу. 
(02 


Выражение (13) описывает сигнал на выходе ключевого элемента, обусловлен- 
ный погрешностью квантования угла Ло, (г). Проанализируем его при отсутствии 


инструментальной погрешности До, что допустимо при [ла . (и) <<1. 
При этом (13) примет вид 


о о и (14) 
ч Г 20 


11 д 
где в, = М,„© — круговая частота ДО... 


Т 
Учитывая, что = (Мо +Ёо)л, где м, = вы ИЕ = [|= | выраже- 
о 
Ч Ч 
ние (14) можно упростить следующим образом: 
[) Т А 
АО, =20 8, »_За(мЕ» )соз1ю, [9 № (15) 
1=1 


где В — период АО. 


Выражение (15) показывает, что погрешность квантования угловой скорости 
на выходе ключевого элемента представляет собой сумму периодической последова- 
тельности импульсов (рис. 4), имеющих высоту 42. и ширину &5Т,, и постоянной 


составляющей —=во@. . Отсюда следует, что среднее значение погрешности кванто- 


вания угловой скорости за период Т, равно нулю. 
АО 


Г 
О-о) 


Рисунок 4 — График зависимости ЛО, (1) - 
Т, — период последовательности импульсов ДО, (+); О, , бо1, — амплитуда 


и ширина импульсов, —6о0, — Постоянная составляющая импульсов 
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На выходе ключевого элемента будут считываться значения ЛО. в дискретные 
моменты времени {= АТ 


ы бо. , А 
во = зодозаи оомь, | 7 (16) 


о 
Поскольку в, =21— и а, = ОГ ‚ выражение (16) можно представить следую- 
4 
щим образом: 


у д 
АО, 7 292 „За (тбо)соз2л7 Борт 


Ч 


(17) 


1 г 
Анализ данного выражения показывает, что для 60 Е | период колебаний 


1 
АО. равен = а для о Е 54 _ ‚ ГДЕ боли =1-бо — Дополнение 6х до 1. 


о И 


1 
Поэтому для &ре | > выражение (17) можно представить в виде 


< 1 
АО и =20 Е. > За (Е )соз2 о, [#-* ь (18) 
1=1 
А, 
где 6, = Е — относительная величина До, выраженная в долях периода. 
обо 
1 
Для ое [5-1 будем иметь 


- 1 
АО; = —20 бои 2, ба (Е лот )соз Ра: |-* + 5 +90 : (19) 


Сопоставление (18) и (19) показывает, что колебания, описываемые (19), инверсны 
колебаниям, описываемым (18), и сдвинуты в сторону отставания на 6, в то время 
как колебания, описываемые (18), сдвинуты в сторону опережения на 00. 


1 ы 
В интервале &5 е и относительная величина периода колебаний ДО, из- 


1 1 
меняется от с до 2, а в интервале 60 =| 0 ‚ что соответствует &се г. з 


относительная величина периода изменяется от 2 до со. 
Переходя к времени / = АТ в (18), будем иметь 


ы т. 
АО и =20 Е, У За (Ес )с03 16.6 | ты) $ (20) 
[= 


2п 
где ©; = = — круговая частота дискретизации отсчетов. 


Отсюда следует, что в случае, когда частота дискретизируемого сигнала пре- 
вышает частоту дискретизации о©„, на выходе получаются отсчеты непрерывного 


№ 
сигнала, имеющего вид исходного сигнала, но в —® +1 раз более низкую частоту. 
[Ф) 
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Поэтому после восстановления запоминанием отсчетов на период будет получен 
непрерывный сигнал АО’ь, имеющий в первом приближении вид исходного сигнала 
ДО, и частоту &о®у. 

При этом импульсы на выходе восстанавливающего элемента имеют примерно 
тот же вид, что и до дискретизации. Отличия обусловлены тем, что в процессе дис- 
кретизации и восстановления осуществляется квантование периода и начального 
смещения исходных колебаний. 

Отсюда следует, что влияние квантования угла при отсутствии инструменталь- 
ной погрешности проявляется в том, что на выходе скоростного канала будем иметь 
сигнал 


®(1)=0+40., (1. 
Поскольку ДО, принимает два значения —Ё00, и (1-Ё,)©,, то ®(!) также 


будет принимать два значения №62, и (№ о По за 


Таким образом, на выходе скоростного канала будут возникать скачки сигнала 
на +6), с частотой = 6 о®, и длительностью 7 (рис. 5). 


©(В д 


№+004 09 


№9 
2п 


бо. 


0 


Рисунок 5 — Зависимость выходного сигнала скоростного канала © (г) : 


0(1) — входной сигнал скоростного канала 


Обозначив © (г) =М, (гро имеем №, (1) = Мо +А№ (1) ‚ где АМ№ (1) — единич- 


ой 


9? 


ные скачки кода длительности Ти частоты в = ы боб а: 


Ей р 
[Ф) 


УФ) 


Чв 


Зависимость АМ№о (г) может быть получена прибавлением к средней со- 
4 


ставляющей 65 


[Ф) 
АМо (1) = 5 ты: 
а 


При этом зависимость ДО, может быть построена по отсчетам ДО, (АТ ), 


описываемым выражением (20), так как 


АО, (1) =АО, (ЕТ) для! | АТ, (к +17 |. (21) 


ав 
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Однако более просто зависимость АМ, (1) может быть получена непосредст- 
венным вычислением №, по формуле 


М, (#)= вы “А — Еп: еетнь, (22) 


и. 


Ч 
Подставив в (22) © =©{ и положив / = АТ , будем иметь 


М, (ЕТ) = М.Е КЕ, +8, |-- Ви (1) +, |, 


7А\ 
где &, =—“. Отсюда находим, что АМ, (АТ) = Ей КЕ, +8, |-- Е (&-1) 80 +, |, 
ча 
а АМ, (1) = АМ» (АГ) для! >| АТ, (к+1)Т ]. 
Согласно расчету, произведенному для случая К = 1,40, бо =0,3 и &, =0,15, 
последовательность АМ (АТ ) имеет вид 
ооо 1 


09010 о отоото ото о отоат 
Отсюда видно, что последовательность состоит из чередующихся периодов 
единичных скачков и содержит 8 периодов длительностью 37 и 4 периода длитель- 


40Т „1 
ностью 47. Средний период равен Г. = ве =3-—Т , что соответствует расчету периода 
ие 3 
12 
по бо. Среднее значение последовательности равно 6: АМ, = 20 = 0,3. 


Таким образом, при отсутствии инструментальной погрешности До преобразова- 
ние постоянной скорости ©» производится с погрешностью квантования Л6<»,, представ- 
ляющей собой последовательность разнополярных импульсов, амплитуды которых равны 


1 1 
ой 800, ДЛЯ Е Е [6 и +60, —(1-2.)0, для ое 54 
Причем длительность импульсов с амплитудой (1 —-0)© „ Равна Т. Их частота 


2п 
в среднем равна &®,, где ®, = = — круговая частота дискретизации. 


Анализ погрешности квантования угловой скорости 
при наличии инструментальной погрешности 


Произведем анализ влияния погрешности ДОу на выходной сигнал скоростного 
канала при наличии погрешности Ло„. Для низких скоростей © погрешность АО, 


описывается выражением (13). Поэтому, заменив в нем №,© на в, ‚ будем иметь 


бо (т, 


АО, (1) =20 80 (1) У`За (мё., (1)) сз, | +С(ЕТ, к (23) 


где о (Г)= 90 (г) = ва | О(1)=О+Ао, (©/). 


Ч 
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Из (23) следует, что ДО, (1) представляет собой последовательность импульсов, 
модулированных по амплитуде, ширине и частоте. При этом амплитуда и ширина поло- 


жительных импульсов будут равны (1-—&, (1))©, и &5 (#7, 


4? 


а частота — ©, ( +Ао,, (01). 


При дискретизации Д<, (1) на выходе ключевого элемента будем иметь 


АО, (ЕТ) =20 60 (кт) У ва (тёо (ЕТ) соз 2 сит)- . (24) 


Учитывая, что выходной сигнал ключевого элемента можно рассматривать как 
результат дискретизации непрерывного сигнала 


ДО, (г) #200 (> (тео (1)) соз/, ты (КТ т Г 


представим (24) в виде 
ло, (1) = ло, (У. 8(1-иТ). 


При этом сигнал на выходе восстанавливающего элемента, осуществляющего 
запоминание отсчетов на период, можно [7] представить как 


АО, (1) =20 Е (>. за, (1))За(п(1Е,+п))х 
хс05 1 тот т 


Анализ данного выражения показывает, что наибольший вклад в АО, (1 ) 


будут вносить составляющие с частотами (1 —п)® а. ыы и И=О для 


ея 
бо (1) 
Ее 5, п=-1 для >< 12, <1. Поэтому ЛО, (г) будет представлять собой последо- 
вательность разнополярных импульсов длительности Т, имеющих амплитуды 
(1-Е. (1))©,‚ 80 (1), и частоту не выше 5 

Зависимость ЛО, (г) при наличии инструментальной погрешности может быть 


построена при помощи выражения (21), где ДО, (АТ ) определяется выражением 


(11), которое с учетом вышеизложенного можно записать в виде 
ло, (Г) =-0,, | ВЕС (АГ) Ег[&о (#1) +6 ((&1)т)] |. — (25) 


Для удобства расчетов представим С(АТ ) как 


С(ЕГ)= в А, эт ("ОКТ+у, )+ к, 


= 


«Штучний 1нтелект» 2013 №3 421 


Габидулин М.А. 


Анализ показывает, что наибольшее влияние на ДО, будет оказывать вторая 
> Ч 
внутришаговая гармоника СКП. Поэтому наибольший интерес представляет оценка 


для ги =2М, , где №, — коэффициент умножения пространственной частоты СКП. 


тах 


Полагая для упрощения расчета 1, =0, будем иметь 


С(кт)= | Аз», зт2 № ОАТ + Ло ]. (26) 
С учетом (26) выражение (25) примет вид 
АО, (АТ) = - | Ре] .К+РЕ[ Аз», т 2, ОЕТ+ А | |- 
—Рг| о (51) +Ег[ Аз, тм, ( ТА ||, (27) 


— АО (ЕТ 
где ло. (АТ) = ый 
4 


Обозначив ф›„ =2М№,ОТ ‚ представим (27) как 
ло, (ЕГ)= | [ЕН Ах, чан, Ао | |- Еи[ о (+1) Лаз, т (#51) фз, Ао | |. 


На рис. 6 приведены графики АО. (АТ } построенные для случаев 6 = 0.3, 


= т 2п 2п 
Д2м, =3, Ао =4, Ф.у 255 и Фу р Из графиков видно, что значения погреш- 
ности ЛО, (ЕТ) не превышают +1. 
АОд ь (9 
А 
1 нЕ 
10т 
25Т „: 
0 5Т 157 
2от/ 
1 
а) 
АО в (0 
Г 
1 г 
|4 
— 
0 5Т 10т 15Т 20Т 
-1 | 
6) 
Рисунок 6 — Графики зависимости ЛО,., (Г) при наличии погрешности До, (1): 
а) для случая &° = 0.3, Дол, =3, Ло =4, ф., =; 
2п 


6) для случая &, = 0.3, Ак, =3, Мо =4, ф.» и 
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Произведенные исследования показывают, что при наличии инструментальной 
погрешности Ла, (1) погрешность ЛО, (АТ) на выходе восстанавливающего элемента 


представляет собой последовательность импульсов, амплитуды которых в зависимо- 
сти от &р, & (АГ) и С(АГ) лежат в пределах +О.. 


Так как ДО та принимает максимальные значения относительно редко, для оп- 


ределения требований к выбору ©, целесообразно оценить среднеквадратическое 
значение погрешности с ( Ао, ) ь 
Произведем оценку с(ло ‚) для случая отсутствия До, (1) . В этом случае по- 
грешность АОчь будет зависеть только от значения & с Е [о, 1) и принимать для каждого 
значения 6, два значения ДО4ы2 С вероятностями р, =6о, р, =1-6о 
ДОцы =1-&5, (28) 
АО» =, (29) 


Из выражений (28), (29) следует, что ДО.» в диапазоне &, равном [0,1), 


может принимать значения, заключенные в интервале (-1, 1) . 
Полагая, что <) распределена в пределах кванта ©, равномерно, определим 
плотность вероятности погрешности квантования ЛО» : 
На участке ДОь Е [0,1) из(28) находим, что 
Ес =1-АОз. (30) 
С учетом (30) можно записать, что 


Р(^О. =| АО», Аб +4ло, || = (Е [+4 ,20|)- р. (31) 


Обозначив плотность вероятности погрешности АО.» через (ло. и выразив 


р, через ДО.» ‚ получим из (31) (ло. =} — ЛО». 
Аналогичным образом можно показать, что для ЛО.» Е [0, —1) 
"(^9. | =1+АОь В 


Отсюда следует, что погрешность АО.» распределена по закону Симпсона и 
представляет собой композицию двух равномерных законов распределения АО и 


АО» . Поэтому среднеквадратическая погрешность Д@,, будет равна 


с(ло.,) = =. (32) 


Таким образом, для достижения требуемой точности на низких скоростях 62 
выбор требований к Да, и <>, должен производиться с учетом выражений (9) и (32) 


О 
Ло, зв у Ч < АО у 
[а [0 /б 


где у — коэффициент, зависящий от значения доверительной вероятности рл 
оценки погрешности. Для р, =0.9 коэффициент у =1.6 [8]. 


. 
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Выводы 


Проанализировано влияние инструментальных погрешностей и погрешности 


квантования на погрешность преобразования угловой скорости и определены требо- 
вания к выбору их основных Параметров цифровых тахометров. построенных на 
основе растровых СКП, и синтезированы точностные и динамические модели, 
позволяющие исследовать их в составе цифровых электроприводов. 
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Рюпа/ Тасйотеегу оп Каяег Кезоег Моуетеи!5 


ТВе агасе апа]у7е4 е еггог$ оЁ тапзКогтайоп ап зупез1$ оЁ ргес11оп @еИа| 
{аспотеег$ оп гафег гезо]уег тоуетет изе4 ю Юпт Фе Ше-зрее4 сваппе! 41еа1 
4нуез оЁ4упапис гапое зрее4 сопёго! ар о 10° [2]. 

То ежепа Ще гапое оЁ уапаЫе зрее4 сопуецегз аге изе4 т Фе Гогт оЁ шсгетеша| 
эгисагез соп$15Н по оЁа гезо[уег, МосК$ е]естоп по@рПег$ ап4 зрайа| йедиаепсу оЁ ве 
ру!5е ип шстетеии [4]. 

ТБе еНесЕ оЁ шзйлитета! еггогз гезо[уег 41$р]асетепе ап4 апошаг тоуетеп оЁ 
Фе даапитаНоп еггог оп Фе ассигасу оЁ фе апошаг уеосИу сопуегзюп соде ш Ч1еиа| 
{аспотеег Базе оп соипНие Фе питбег оЁ шагу14иа| шсгетеп ог апешаг тоуетепе 
ОР Фе соде Юг а Нхе4 репо4 оЁ Ите ап4 ассигасу гедитетеп 4ейпе4 Бу Ше сво1се оЁ 
(5$ рагатебегв. 

Г 15 езаБИ$Бе4 фа Фе шзбглитена! еггог ог е тасвотеег 1$ затр!е4 а{ е ите 
{фе депуайуе шзубитепе еггог$ гезо!уег ФаЁ Ваз раззе гоиеН Ше ауегасше ап 
тедисте Фе е]етеп Фа 4езсте Фе апзе оЁ Ч@етепнайоп ул гезресе © Ите апа 
тетогу затр/ез рег сусе а{ %ер арргохипайоп. 

Зупе$17е ассигасу ап4 Чупапис то4е|$ №0 ехроге е Чоиа| ‘асвотаег$ оп 
газ{ег гезо|уег тоуетеп$ ш 41еИа| апуез. 
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